Violation of area- law scaling for the entanglement entropy in spin 1/2 chains 
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Entanglement entropy obeys area law scaling for typical physical quantum systems. This may 
naively be argued to follow from locality of interactions. We show that this is not the case by con- 
structing an explicit simple spin chain Hamiltonian with nearest neighbor interactions that presents 
an entanglement volume scaling law. This non-translational model is contrived to have couplings 
that force the accumulation of singlet bonds across the half chain. Our result is complementary to 
the known relation between non-translational invariant, nearest neighbor interacting Hamiltonians 
and QMA complete problems. 
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I. INTRODUCTION 

Ground states of relevant physical Hamiltonians carry 
quantum correlations which decrease with distance. For 
instance, a two-point correlation function is expected to 
fall off exponentially with the separation of points in the 
presence of a mass gap or algebraically at critical points. 
This amount of quantum correlations is in correspon- 
dence with an area law scaling of entanglement entropy. 
To be precise, the entanglement entropy, defined by 



S{pa) = -tr {pA^ogpA) 



(1) 



where pA is the density matrix of the region of space A 
considered, scales as the boundary of A. 

In contradistinction to the above situation, random 
quantum states are known to carry volume law entropy 
[l||. Therefore, typical Hamiltonians produce grounds 
states which are not generic. Indeed, relevant physics 
corresponds to a small corner of the total Hilbert space 
of a quantum system. This observation is crucial to un- 
derstand recent efforts to simulate quantum states with 
tensor networks Such approximations are able to 

accommodate area law scaling for the entropy. 

It is then important to understand precisely what are 
the properties that a Hamiltonian must obey so as to 
produce a ground state which only displays area law en- 
tanglement. A first heuristic approach suggests that en- 
tanglement decreases at large distances because interac- 
tions are local. That is, a local degree of freedom inter- 
acts with its neighbor and gets entangled with it. This 
second degree of freedom interacts in turn with a fur- 
ther one. This sequence of interactions would eventually 
entangle far separated degrees of freedom, though the 
strength of interactions would only manage to get the 
standard correlations we find in Nature. On the other 
hand, it is unclear whether interactions could be con- 
trived to achieved larger entangled states. The role of 
translational invariance is then critical. 

Some results for one-dimensional systems are well- 
established. In one dimension, if a system obeys local 
interactions and it is gapped, area-law always emerges 
0. On the other hand, if the system is at a critical 



point, and therefore gapless, a logarithmic divergence is 
encountered. This logarithmic scaling of the entangle- 
ment entropy is explained by conformal field theory |lOl - 

H. 

In Refs. (TBI . [TgI . infinite translational invariant 
Fermionic systems of any spatial dimension with arbi- 
trary interactions are considered. For such systems, it is 
shown that the entropy of a finite region typically scales 
with the area of the boundary times a logarithmic cor- 
rection. 

Although there has recently been further progress on 
this topic in higher dimensions [iTl - ITgt , the necessary and 
sufHcieirt conditions for an area-law have not been defined 
yet. An explicit example of a system where area-law is vi- 
olated is presented in Ref. [23|. It is shown that a one di- 
mensional non-translational invariant system composed 
of 12-level local quantum particles with nearest neighbor 
interactions presents a ground state that carries a volume 
law scaling of entanglement. In particular, it is proven 
that the problem of approximating the ground state en- 
ergy of such system is QMA-complete. This precise ex- 
ample shows that a quantum computer could not sim- 
ulate any one dimensional system, and, moreover, that 
there exist one-dimensional systems which take an expo- 
nential time to relax to their ground states at any temper- 
ature, making them candidates for being one-dimensional 
spin glasses. 

The issue addressed in this work is to study how sim- 
ple can be a quantum system to give a highly entan- 
gled ground state. In particular, we show that a simple 
spin 1/2 model with nearest neighbors interactions with 
a suitable fine tunning of its coupling constants can have 
a ground state with a volume law scaling for the entangle- 
ment entropy. Our proposal is based on the translational 
symmetry breaking, hence, this makes that the area-law 
violation can not be maintained for any bipartition of the 
system. Nevertheless, it will be shown that the average 
of the entanglement entropy over all the possible posi- 
tions of the block fulfills a volume-law. Our results are 
presented in the following way. We first review the real 
space renormalization group technique which brings the 
fundamental intuition on how to build an area law vio- 
lating Hamiltonian. We then turn to solve the proposed 
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Hamiltonian using its exact diagonalization. where the 
final step can be taken both in perturbation theory or 
numerically. We also illustrate the real space renormal- 
ization idea in an Appendix. 



II. REAL SPACE RENORMALIZATION GROUP 

A. Introduction to real space Renormalization 
Group approach 




Figure 1: Diagram of a random singlet phase (a) and the 
concentric singlet phase (b). Each spin forms a singlet pair 
with another spin indicated by the bond lines. 



Real-space Renormalization Group (RG) approach was 
introduced in Ref. [2l| generalizing the works presented 
in Ref. . It is a method suited for finding the effective 
low energy Hamiltonian and the ground state of random 
spin chains. The couplings have to satisfy the hypothesis 
of strong disorder, i. e. the logarithm of its probability 
distribution is wide. Under such conditions, the ground 
state of the system can be very well approximated by 
a product state of singlets whose spins are arbitrarily 
distant. 

Let us review the real-space RG method for the inho- 
mogeneous XX model case 

1 ^ 

fc=2E'^'K^m+'^M+i) • (2) 

i=l 

First, we find the strongest bond Ji ^ Ji+i,Ji-i and 
diagonalize it independently of the rest of the chain. Ac- 
cording to the previous Hamiltonian, this leads to a sin- 
glet between spins i and i + 1 (see appendix |X| . There- 
fore, the ground state at zeroth order in perturbation 
theory respect the couplings Ji-i and Ji-^-i is 

|V'^°^> = lV'.<.)l^->lV'.>.+i> (3) 

where ji/j-) = (|01)i_i+i — |10)i,i+i) is a singlet state 

between the spins i and i + I, and \4'x<i) and \4'x>i) 
correspond to the states of the rest of the system. 

In order to compute corrections to the ground state at 
higher orders, we use perturbation theory as it is shown 
in appendix [Xj This leads to an effective interaction be- 
tween the distant spins j — 1 and i + 2 with an effective 
coupling 

i.-M+. - '-^^ ■ (4) 

In summary, real space RG integrates out two spins and 
reduces the Hamiltonian energy scale. Notice that this 
new effective low energy Hamiltonian couples the spins 
i — \ and z -|- 2, therefore, it has non-local interactions as 
seen from the original Hamiltonian. Iterating this proce- 
dure for a XX model with random couplings, it is seen 
that the ground state can be described by a random sin- 
glet phase, i. e. each spin forms a singlet pair with an- 
other one (see Fig. [T^). Most pairs involve nearby spins, 
but some of them produce long distance correlations. 



In Ref. [l^, real-space RG was used to show that, for 
random spin chains where the ground state is a random 
singlet phase, the entanglement entropy scales logarith- 
mically at the critical point as in the homogeneous case. 
That is, 

SL^^log^L, (5) 

where c = c In 2 is an effective central charge propor- 
tional the central charge for the same model but without 
disorder c. This ana.lytical result has later been check 
numerically in Refs. [mHII. 



B. Area-law violation for the entanglement entropy 

Let us now tune the couplings Ji of our XX model in 
such a way that the entanglement entropy of the ground 
state of the system scales with the volume of the block 
of spins. An easy way of achieving this is to generate a 
ground state with a concentric singlet phase as it is shown 
in Fig. [T|3. We see that the system is in a product state 
of distant singlets between the positions N/2 — (z — 1) 
and A^/2 + i for 1 < i < A^/2. It is trivial to see that the 
entanglement entropy of this configuration would scale 
with the size of the block, since it merely corresponds to 
the number of bonds cut by the bipartition (see Fig. [5^) . 

Let us note that the entanglement entropy for con- 
centric blocks would be as it is shown in Fig. ^p. As 
translational invariance of the system is broken, the en- 
tanglement entropy of a block not only depends on the 
size of it but also in its position. 

In order to measure how entangled is a state for non- 
translationally invariant systems, it is useful to introduce 
the average entanglement entropy over all the possible 
positions of the block, that is 

N-L 

- ]^ E SL{^) (6) 

i=l 

where S'i(i) is the entanglement entropy of the block of 
size L from the i-th spin to the {i + L)-th one. 

According to the previous definition, the average en- 
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J^ 




(b) 




Figure 2: Diagram of the entanglement entropy scaling for 
the concentric singlet phase. The entanglement entropy grows 
maximally if we take blocks at one extreme (a) and is zero 
if the blocks are centered at the middle of the chain. This is 
an explicit example that in the non-translationally invariant 
systems the entanglement entropy depends on the position of 
the block. 



tanglement entropy of the concentric singlet phase reads 



Sr 



2{N - L] 



L. 



(7) 



Although for the concentric singlet phase the average 
entropy losses its linear behavior for large blocks, L ^ 

^1 — N, it always fulfills the condition Sl > \L- 

Thus, the concentric singlet phase represents a simple 
and explicit example of area-law violation of scaling of 
the entanglement. 

The aim of our work is to tune the coupling constants 
Ji of the XX model, such that, the concentric singlet 
phase becomes the ground state of the system, and, in 
this way, to obtain an explicit example of a Hamiltonian 
with nearest neighbor interactions of spins that violate 
the area-law scaling of entanglement. 

Due to the symmetry of the state that we pretend to 
generate, let us consider a XX chain of N spins where 
the central coupling between spins N /2 and N/2 4- 1 is 
Jo and the rest of them are chosen as follows 



J. 



N/2+i,N/2+i+l 



J. 



JV/2-i,JV/2-i+l 



(8) 



where 1 < i < ^ — 1 and the coupling JN/2±i connects 
the spins N/2 ± i and N/2 ±i + l. 

We are going to use real space renormalization group 
ideas in order to see at which values we have to tune 
the coupling constants, such that, the concentric sin- 
glet phase becomes the ground state of the system. If 
Jq ^ Ji, in the low-energy limit, an effective interaction 
between the spins N/2—1 and N/2 + 2 appears. We label 
this effective coupling as Ji and, according to Eq. Q, it 



2Ji- 



Figure 3: Diagram of the formation of an effective coupling 
Ji if the condition Ji-i 3> Ji is fulfilled. 



reads 



Ji 



2Jo 



(9) 



Then, if Ji 3> J2, the effective low-energy Hamiltonian 
will have an effective bond between the spins N/2~-2 and 
N/2 + 3. We would like to proceed in this way in order 
to generate iteratively the concentric singlet phase. 

Thus, if the condition Ji ^ Ji+i is fulfilled in general, 
where Ji is defined by 



J. 



J 



2J,- 



(10) 



we expect that the ground state of the system is the 
concentric singlet phase. 

Specifically, if we impose that Ji = eJi-i for any i, 
such that it is always possible to apply Eq. ((TO)) . we see 
that the couplings Ji must decay very rapidly 



• 



(11) 



In general, we are going to study chains with couplings 
that decay 



(12) 



where a{i) is a function that is monotonically increasing. 
If a(i) ~ i^, we would have a Gaussian decaying. 

Next; we want to solve the XX model with the coupling 
constants defined in Eq. ((T^), and study how the entan- 
glement entropy scales depending on the kind of decay 
law for the couplings. 



III. SOLUTION OF A SPIN MODEL AND ITS 
ENTANGLEMENT ENTROPY 

Let us consider a finite spin chain with nearest neigh- 
bor couplings Jf , Jf and an arbitrary transverse mag- 
netic field Xi in each spin. This system is described by 
the Hamiltonian: 



H 



1 ^ 

4i: 



' TX ^X 
,Ji CTj CTj + l 



1=1 



N 

E- 

i=l 



(13) 



where L is the size of the system and cr^'^ are Pauli- 
matrices at site i. The XX model presented before is a 
particular case of this Hamiltonian (|13p for Jf = J^ and 
A, = Vi. 
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A. Jordan- Wigner transformation 

The essential technique in the solution of H is the 
well-known mapping to spinless fermions by means of 
the Jordan- Wigner transformation. First, we express the 
spin operators cr^'^'^ in terms of fermion creation (anni- 



hilation) operators cj (ci): c| = a,^ exp 



3 '^jS 



and Ci ~ cxp 



Ei-l + - 



a, , where a - 



ia^)/2. Doing this, H can be rewritten in a quadratic 



form in fermion operators 

TV 1 ^ 

2 



ij"=l i,j = l 

where the matrices A and i? are defined by 



(14) 



A, 



2XA,j + {Jf + Jf )<5,+ij + ( j; + Jf )5.,,+i 



(15) 



with l<i,j <N. 

B. Bogoliubov transformation 



In a second step, the Hamiltonian is diagonalized using 
a Bogoliubov transformation 



1=1 



(16) 



where the $fc and are real and normalized vectors: 



E: '^'(O = E: = 1, so that we have 



N 



(17) 



fc=i 



The fermionic excitation energies, A^, and the compo- 
nents of the vectors, and 'i'k, are obtained from the 
solution of the following equations: 



{A - B)^k = A,.*fc 
{A + B)*fc = Afc$fc . 



(18) 
(19) 



It is easy to transform them into an eigenvalue problem. 



(A + B)(A-B)$fc = A2$fc 
{A~B){A + B)^k^Al-^i,, 

from where A^, $fc and 5'^ can be determined. 

C. Ground State 



(20) 
(21) 



allows us to restrict ourselves to the sector corresponding 
to Afc > 0, fc = 1, 2, . . . , Af. Thus, considering Eq. ^ 
and the fact that all Afc are positive, the ground state is 
a state \GS) which verifies. 



(22) 



In practice, what we do to restrict ourselves to the 
sector of positive Ak is to determine $fe and Ak by solving 
Eq. (HOI), and calculate ^-fc = ^(A - B)^k- 



D. Computation of Von Neumann entropy 

Following Refs. [Til. [27l. [28| . the reduced density matrix 
PL = tr n^l\GS){GS\ of the ground state of a block of 
L sites in a system of free fermions can be written as 



PL = Ke 



(23) 



where k is a normalization constant and % a free fermion 
Hamiltonian. 

Let us very briefly justify why the density matrix must 
have this structure. First, notice that the Hamiltonian 
defined by Eq. ([T4)) has Slater determinants as eigen- 
states. Thus, according to Wick theorem, any correlation 
function of the ground state (or any other eigenstate) can 
be expressed in terms of correlators of couples of creation 
and annihilation operators. For instance, 

(4cLcfcC/) = {cici){c\^Ck) - (4cfe)(cJ„Q) + {cicl,){ckci) . 

(24) 

If all these indices belong to a subsystem of L sites, the 
reduced density matrix must reproduce the expecta- 
tion values of the correlation functions, i. e. 

tr (pLclcl^CkCi) = tr (plcIci) tr (picJ^Cfc) 

-tr {pLclck)tT {pl4„,ci) (25) 
+ tr {pl ) tr (pLCkCi) ■ 

This is only possible if p^ is the exponential of an op- 
erator Ji which also contains creation and annihilation 
processes, i. e. 



n 



L 

yi ^Mji-i>-j . 2 



h.c. 



(26) 



We can diagonalize this Hamiltonian H by means of 
another Bogoliubov transformation 

(27) 

where Vk{i) and Ufc(i) are real and normalized. Then, the 
Hamiltonian reads 



In Eqs. ([T^ and ([T5)) . we realize that transforming ^k 
into — $fc (or vj/j, into — ^fc), Ak is changed to — A^-. This 



n = J2ek^Uk 

k=l 



(28) 
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where and are the creation and annihilation oper- 
ators of some fermionic modes. In terms of these modes, 
the density matrix is uncorrelated and can simply be 
expressed as 



where 

Pk = 



1 + e-"" 



PL 



1 



2 " 




(29) 



(30) 



In the previous equation, the new parameters Vk have 
been introduced in order to ensure the normalization of 
Pk, tr (yOfc) = 1. This way of expressing pk will be useful 
next. 

Thus, the entanglement entropy of the density matrix 
Pi is merely the sum of binary entropies 



5(i) = ^5(p,) = ^i/ 



fc=i 



k=l 



l + Vk 



(31) 



where -ff(p) = — plog2p— (1 — p) log2(l —p) is the binary 
Shannon entropy. 

In order to determine the spectrum of p^, let us con- 
sider the correlation matrix, 



(G5|(ct -c„)(4+c„)|G5). (32) 



Notice that the matrix G can be computed using the $fc 
and vectors. 



N 



Gm.^n = - ^*fc(m)$fc(n). 



(33) 



k=l 

.J, 



where the correlations {rikVq) — ^kq and {rjk'>]q) — have 
been considered. 

In the subspace of L spins, G is completely determined 
by the reduced density matrix. To avoid any confusion, 
let us define T = G(l : L, 1 : L) as the L x L upper-left 
sub-matrix of the correlation matrix G. Then, T can be 
expressed in terms of the expected values {^l£,q), 

L 

k,q=l 
L 

= ^Ufc(i)wfc(j>fe , (34) 

k=l 

where the i and j indices run from 1 to L. This equation 
leads to the relations, 

TUg = l^gVg (35) 

T'^Vq = VqUq , (36) 

that can be translated to the eigenvalue problem 

T^TUq = VgUq (37) 

TT^Vq = ulvg . (38) 

Once the Vq variables are computed, we can determine 
the entanglement entropy by means of Eq. pip . 



E. Summary of the calculation 

To sum up, let us enumerate the steps that we have to 
follow in order to calculate the entanglement entropy of 
a block L. 

1. Write down the matrices A and B in terms of the 
couplings of the Hamiltonian according to Eqs. 

(Ha. 



2. Determine A^, $fe and '^k by solving the eigenvalue 
problem from Eq. (|20p . 



3. Calculate the correlation matrix G defined in Eq. 
(EH). 

4. Take the sub-matrix T and to determine the eigen- 
values Uk from Eq. 



5. Compute the entanglement entropy by means of 
Eq. (EJ). 



IV. EXPANSION OF THE ENTANGLEMENT 
ENTROPY 

We would like to tune the coupling constants of the 
Hamiltonian p^ . such that the scaling of the entangle- 
ment entropy of its ground state violates the area law. 
The entanglement entropy only depends on the variables 
Vk- Then, we can separate the Shannon entropy of the 
probabilities ^^p- into 



H 



l + Vk 



1 - h{uk) 



(39) 



where h{x) = — ilog(l — x^) — ^ (l+f)' ^ positive 
function. Thus, the entanglement entropy reads 



S{L) = L -Y,K^k) ■ 



(40) 



fe=i 



Notice that the scaling of the entropy only depends on 
the sum X]fc=i ^{^k)- More concretely, we can define the 
parameter. 



1 ^ 

L— ^OO Li ^ — ^ 



(41) 



k=l 



that describes the asymptotic behavior of the scaling of 
the entropy for large blocks: 



• /3 = 0: maximal entanglement, 

• /3 < 1: volume-law, 

• B = 1: sub- volume-law. 
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Let us focus on the case ^ ~ 0. Let us analyze if it 
is possible to design a spin chain with nearest neighbor 
interactions whose ground state is maximally entangled. 
First, we realize that /3 is strictly zero if and only if all the 
variables Vk = 0- Thus, if we want to consider small de- 
viations of the maximally entangled case, we can assume 
that Vk ~ and expand /3 in series of 



1 



(42) 



k=l 



Considering Eq. psp . we can express (3 in terms of the 
matrix-elements of T, 



/? = hm tr (TT^) = lim V - . 



(43) 



Let us notice that to fulfill this condition requires that 
the average of the matrix-elements of T tend to zero for 
large L, 



1 



L 

lim — y 



0. 



(44) 



If we assume a smooth behavior for the matrix-elements 
of T, according to Eq. ((43|) . they must decay faster than 
the inverse square root function. 



1 



(45) 



such that /3 = 0. 

In conclusion, in order that the entanglement entropy 
scales close to the maximal way, the matrix-elements of 
T matrix have to be very close (or decay rapidly) to zero. 
If this is the case, the entanglement entropy can be sim- 
plified to 



S{L)=L-\\T\ 



F I 



(46) 



where ||T||i? is the Frobenius norm of T, defined by 
\\T\\f = {T^T). _ ^ 

Let us now study if it is possible to tune the coupling 
constants of a spin chain in order that ||r||_F is strictly 
(or close to) zero. The possibility of having a null T 
is discarded because it cannot be achieved with nearest 
neighbor interactions models. Despite this, there is a 
wide freedom to tune the coupling constants such that 
the matrix-elements of T fulfill condition (|43|) . This ar- 
bitrariness makes very difficult to specify the shape of the 
distribution of coupling constants in order that area-law 
is violated. With this aim, we can exploit the idea of real 
space Renormalization Group presented before. 



NUMERICAL RESULTS 



We can follow the steps described in Sec. lIIIEl in order 
to calculate the entanglement entropy of the XX chain 



12 

10 - 



4 
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(a) 



entropy 
trace-norm 
first-term 
0.697367*x+0 .309923 
,097*x'\1 .32)1+2.1 - 
i9'f1(x)-0.9n2fxl - 



(b) 



entropy 
trace-norm 
first-term 
0.9805*x-0.34 
0.841*x-1.18 ' 
1.0'f1(xl-1.0'f2lxl - 



Figure 4: Scaling of the entanglement entropy of a block of 
contiguous spins with respect to the size of the block L for the 
ground state of a XX model with couplings that decay: (a) in 

2 

a Gaussian, J„ — e~" , and (b) in an exponential, J„ = e~", 
way. The magnetic field is set to zero. The Frobenius norm of 
the T matrix related to this system and the sum S{L) + 1 |r| |f 
are also plotted. This allows us to check how accurate is the 
approximation (|46p both in the case of Gaussian decay (a) 
and in the exponential one (b). 



presented in Sec. Ill Bl and check if the entanglement en- 
tropy grows linearly with the size of the block. 

This XX model is characterized by having the strongest 
bond in the middle of the chain, Jq, while the value of 
the rest of bonds J„ decrease rapidly with the distance 
n to the central one. In particular, we have studied two 
different kinds of decay for the coupling constants J„: (i) 
Gaussian decay, Jn = e~" , and (ii) exponential decay 
Jn = e^". Let us notice that due to the rapidly decaying 
of the coupling constants and the finite precision of the 
computer, we can only consider small systems. 

In Fig. Ufa), the entanglement entropy is plotted for 
the Gaussian case. As we expected scales linearly with 
the size of the block L with a slope practically equal to 
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one. Notice that although the slope is 1 for large blocks, 
the entanglenient is not the maximal due to the non- 
linear behavior of the entropy for the smaller ones. This 
can be better understood analyzing if the approximation 
of the previous section given by Eq. (pS)) is fulfilled. With 
this aim, the square Frobenius norm HTHln and the sum 
S'(L)-|-| |T|||, are also plotted. In fact, we observe that the 
sum S{L) + ||T|||n coincides with the maximal entropy, 
as Eq. suggests. 

The same plot can be realized for an exponential decay 
of the coupling constants, see Fig. IH^b). In this case, 
although the entanglement entropy also scales linearly, 
its slope is less than one. Thus, we observe a volume 
law, but the entropy is not maximal. Therefore, Eq. (H^ 
is not fulfilled in this case. We can, actually, see that the 
Frobenius norm ||r||j;^ increases linearly with L instead 
of saturating to a small value. 

We have repeated the same computations for the same 
kind of decays but other basis. The same behaviors for 
the Gaussian and the exponential cases have been ob- 
tained. For the Gaussian faster decay implies 
a saturation to a smaller value for that is, a 
closer situation to the maximal entropy. For the expo- 
nential case, ||T||i? continues increasing linearly but with 
a smaller slope. 



VI. CONCLUSIONS 

We have constructed a one dimensional system com- 
posed by spin-i particles with nearest neighbor interac- 
tions with a entanglement entropy of the ground state 
that scales with the volume of the size of the block. 

This results further confirms that violations of area law 
scaling for entanglement entropy are possible for local 
interacting Hamiltonians. Furthermore, such a behavior 
is found possible for spin 1/2 degrees of freedom. The 
price to be paid for violating area law scaling is to break 
the translational symmetry of the system. Indeed, in 
Refs. [2^ [sOl, it is shown that, although translation- 
invariant one dimensional states give rise to arbitrary fast 
sublinear entropy growth, they cannot support a linear 
scaling. 

Let us also recall that two-local Hamiltonian problems 
have been proven to provide QMA-complete problems. 
To be more precise, the problem of finding out whether 
the ground state energy of a two-local Hamiltonian is 
larger than a or smaller than b, where |a — 6| > 0(l/n) is 
QMA-complete. We may further argue that an eflicient 
classical simulation of such a problem is likely to be im- 
possible, otherwise we could solve any NP-complete prob- 
lem by just simulating Quantum Mechanics on a classical 
computer. The obstruction to obtain faithful simulations 
of quantum mechanical systems is in turn related to the 
amount of entanglement found in the system. Thus, ex- 
ponentially large entanglement should be found in some 
one-dimensional quantum systems. Our results is some- 
how completing this idea. Even spin 1/2 chains can pro- 



duce highly entangled states if couplings are adequately 
timed. 



Appendix A: Real space renormalization group in a 
XX model of 4 spins 

We consider first a simple XX model with only 4 spins 
and couplings {A, 1, A}. We can rewrite the Hamiltonian 
of the system as a perturbation theory problem, 

H = Ho + XV, (Al) 

where, 

Ho^a^a-^ +alal , (A2) 

and 

V = (7^ (7^ + aY + (7^ + (tJ . (A3) 
The eigenstates of Hq are 

|V^+) = -^ (101)23 +110)23) 
1^0) = 100)23 

iV'i) = 111)23 (A4) 
IV-) = ^(101)23-110)23) 

with eigenvalues +2, 0, and -2 respectively. We are in- 
terested in study what happens to the ground state (GS) 
of the Hamiltonian H when the perturbation XV is intro- 
duced. The ground state of Ho is degenerate and form 
a subspace of dimension 4. In particular, we choose the 
set of vectors {|m)} = {|0)i|i^_)23|0)4, |0)i|V'-)23|l)4, 
|l)i|V'-)23|0)4, |l)i|'0-)23|l)4} as a basis. 

We expect that the perturbation removes the degener- 
acy in the sense that there will be 4 perturbed eigenstates 
all with different energies. Let us call them {|/)}. As A 
goes to zero, \l) tend to 1^'°-') which are eigenstates of Hq, 
but which in general will not coincide with | m) . 

According to perturbation theory, let us expand the 
eigenstates and the eigenvalues of H in powers of A, 

10 = |;W) + Al^i)) + A^lO^)) + 0(A3) (A5) 

and 

El = E^°l + \e[^^ + \^e\^^ + 0{X^) . (A6) 

Notice that the zero order term in the energy expan- 
sion does not depend on Z, since the ground state of the 
non-perturbed Hamiltonian is degenerate. Substituting 
the previous expansions into the Schrodinger equation, 
[Hq -\- W) \l) = Ei\l), and equating the coefficient of var- 
ious powers of A, we obtain a set of equations that will 
allow us to find the corrections to the perturbed eigen- 
states and eigenvalues. 
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At zero order in A we recover the trivial non-perturbed 
Schrodinger equation. If we collect terms of order A, we 
get 



(A7) 



In order to calculate the first correction to the energy, 
we project the previous equation (|A7p to the degenerate 
ground state subspace 



where al^ are the coefficients of jZ^*^^) = X^m'-'^ml"^) 
pressed in terms of the basis |m). Notice that Eq. (jA12p 
is a diagonalization problem. For our particular Hamil- 
tonian, it takes the form 



/I 0^ 

110 

110 

VO 1/ 







/ai 








"3 1 










\a4 



(A13) 



(A8) 



Now the degeneration is completely broken and the per- 
turbed ground state becomes 



where Kn.m' = (m|l/|m') is the projection of the interac- 
tion to this subspace. In our particular case, the matrix- 



elements Vm 



for all TO and to', hence, 



(1) 



V 



I . This means that the degeneration is not broken at first 
order in A and forces us to consider the second order. 

Ho) = {V- eI'^) 1^(1)) - i?P>|Z(°)) . (A9) 

We proceed as previously and project this equation to 
the degenerate ground state subspace, 



{m\V-El'^\l^^^)=EP{m\l^°^) 



(AlO) 



\GS) = |V-'-)i4|V'-)23 - A^(llOOl) + 10110)) + OiX"), 

(A14) 

with 



Egs ^ -2 + 0{\^) . (A15) 
We can obtain an effective Hamiltonian by projecting the 
original one into the subspace of lower energy formed by 



\2 

Heff = PHP^ = -2+— (2 + a^af + aj aj) , (A16) 



From equation (|A7P we can compute the first order cor- 
rection to the eigenstates \l), 



(fc(o)|y|;(o)) 

I' / ^(0) j^(0) 



k^GS ^GS 



(All) 



where |fc(°)) are the Ha eigenstates that do not belong to 
GS. Now, we substitute this into (jA9p and get an equa- 
tion for the second order correction to the energies and 
the states |Z(°') 



E 



TO|y|fcW)(fcW|F|TO) , (2) , 



p(0) _ p(0) 
m',k ^GS ^k 



(A12) 



where P = Y.m l^^"')^^"'! and e {|V-_)i4|^-)23, 

|'/'o)l4|V'-)23, |V'l)l4|?/'-)23, |V'+)l4|V'-)23}- 
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